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The direction of the steady-state heat currents across a generic quantum system connected to multiple baths 
may be engineered so as to realize virtually any thermodynamic cycle. In spite of their versatility such contin¬ 
uous energy-conversion systems are generally unable to operate at maximum efficiency due to non-negligible 
sources of irreversible entropy production. In this paper we introduce a minimal model of irreversible absorp¬ 
tion chiller. We identify and characterize the different mechanisms responsible for its irreversibility, namely 
heat leaks and internal dissipation, and gauge their relative impact in the overall cooling performance. We 
also propose reservoir engineering techniques to minimize these detrimental effects. Finally, by looking into a 
known three-qubit embodiment of the absorption cooling cycle, we illustrate how our simple model may help 
to pinpoint the different sources of irreversibility naturally arising in more complex practical heat devices. 

PACS numbers: 05.70.-a,05.30.-d,07.20.Pe 


I. INTRODUCTION 

Almost two centuries ago, Sadi Carnot established the max¬ 
imum efficiency attainable by heat engines in his ground¬ 
breaking work on the motive power of heat [1]. He did so 
by considering a simple reversible model (the Carnot cycle), 
which was general enough to account for any specihc realiza¬ 
tion of a steam engine in the limit of quasistatic operation. 
Realistic devices, however, need to be operated in hnite time 
and are, therefore, intrinsically irreversible. Thus, in order 
to address the practical problem of performance optimization 
of irreversible engines (or refrigerators), better benchmarks, 
tighter than the Carnot bound, must be obtained. This is 
the central question of finite-time thermodynamics [2], and 
addressing it requires to introduce suitable models, simple 
enough to remain analytically tractable, and still contain all 
the relevant physics of an itTeversible device. 

The hrst attempts relied on endoreversible models, which 
assume that hnite-rate heat transfer effects yield the leading 
contribution to the total irreversible entropy production, and 
disregard completely other naturally occurring irreversible 
processes such as internal dissipation, friction or heat leaks. 
This approximation, originally intended to account for sim¬ 
ple turbine models [3-5], has been generically used in the de¬ 
scription of heat engines and refrigerators [6-8]. However, a 
careful analysis shows that the operation of commercial de¬ 
vices is indeed dominated instead by internal dissipation and, 
to a lesser extent, heat leaks, which advises against its use [9]. 

Finite-time thermodynamics has also been applied to quan¬ 
tum energy conversion cycles, i.e. those embodied in individ¬ 
ual quantum systems with discrete energy spectrum [10-12]. 
These nanoscale heat devices were hrst studied back in the 
late 1950s [13-15], when it was pointed out that a three-level 
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maser could be thought-of as the smallest realization of a heat 
engine, or a compression refrigerator. Interestingly, a weakly 
driven three-level maser is, in good approximation, endore¬ 
versible [16], while a three-level absorption chiller or heat 
transformer (i.e. driven by heat rather than work) [17, 18] 
is strictly so. The optimization of this model has been thor¬ 
oughly studied so as to establish universal endoreversible per¬ 
formance bounds from a full microscopic description of the 
system-baths interactions [19-23]. 

Very recently, several experiments have been proposed to 
realize quantum energy conversion devices, such as a heat en¬ 
gine made up of a single trapped ion [24, 25], laser-driven 
two-level chillers [26-28]; or even absorption chillers real¬ 
ized in optomechanical setups [29], arrays of coupled quan¬ 
tum dots [30], or super-conducting qubits [31]. Hence, the 
interest in nanoscale heat devices is starting to transcend the 
theoretical understanding of the foundations of thermodynam¬ 
ics as their hrst practical applications start to be envisaged. 

Just like in the macroscopic realm, realistic nanometer¬ 
sized engines and refrigerators are generally not strictly en¬ 
doreversible [11]. Therefore, in order to predict characteristic 
parameters such as the maximum cooling coefficient of perfor¬ 
mance (COP), or to hnd generic design prescriptions compat¬ 
ible with operation at maximum power (or cooling rate), one 
needs to come up with a prototype model, complex enough to 
encompass all the relevant lossy mechanisms other than the 
unavoidable hnite-rate heat transfer effects. 

For instance, quantum Otto cycles performed infinite time 
on interacting two-level atoms were proposed for the study of 
‘friction’ [32-35], which was identihed with the irreversible 
degradation of residual quantum coherence. When it comes 
to continuous models (i.e. those operating under steady-state 
conditions), it is known that e.g. a weakly-driven three-level 
chiller does suffer from ‘heat leaks’ once the parasitic cou¬ 
pling between its driving transition and the environment is 
accounted for [17]. As shown by Kosloff and Levy in [10], 
even perfectly isolated three-level chillers deviate from en- 
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doreversibility under strong driving. This is due to an essen¬ 
tially different type of losses, which result from an internal 
competition between alternative decay paths. This losses are 
accompanied by heat leaks. Interestingly, both internal losses 
and heat leaks are also present in absence of external driv¬ 
ing or imperfect thermal insulation. Indeed, multi-terminal 
autonomous devices, i.e. those consisting of mutually inter¬ 
acting contact ports locally addressed by independent heat 
baths, operate irreversibly as a result of ‘delocalized dissipa¬ 
tion’ [20]. 

In spite of the efforts devoted to the study of irreversibil¬ 
ity in nanoscale continuous devices, a complete understand¬ 
ing of its basic underlying mechanisms and their interplay 
is still missing. In this paper we intend to hll that gap by 
showing how a generic quantum heat engine or refrigerator 
can be thought-of as a ‘multi-effect’ device [22], made up 
of detuned stages that may work cooperatively and, crucially, 
also compete with each other [10]. We genetically term this 
competition ‘internal dissipation’. Heat leaks appear then as 
a mere by-product, since the combination of such stages fa¬ 
cilitates parasitic energy flows between heat baths. Addition¬ 
ally, imperfect insulation at the thermal contacts introduces 
undesired system-bath couplings that result in a direct bath- 
to-bath heat transfer or ‘thermal short-circuit’. To carry out 
our analysis, we propose a minimal irreversible absorption re¬ 
frigerator which reduces to a double-stage absorption chiller 
[18], with detuned three-level stages. We further quantify and 
compare the individual contributions of internal dissipation 
and heat leaks to the total irreversible entropy production, and 
even propose reservoir engineering techniques to partly sup¬ 
press these detrimental effects and thus, increase the cooling 
performance. Finally, we show how the basic lossy mecha¬ 
nisms present in our simple model also show up naturally in 
more complex practical multi-terminal interacting absorption 
chillers [20]. 

This paper is structured as follows: In Sec. II we intro¬ 
duce our minimal irreversible model. Its steady-state solu¬ 
tion is discussed in detail in Sec. Ill: In particular, we propose 
a meaningful breakdown of its steady-state heat currents in 
Secs. Ill A and III B. We also discuss the phenomenon of ther¬ 
mal short-circuit in Sec. Ill C, and propose strategies to partly 
suppress internal dissipation and heat leaks in Sec. HID. The 
diagnosis of the more complex three-qubit interacting model 
is tackled Sec. IV. Finally, in Sec. V we summarize and draw 
our conclusions. 


II. THE FOUR-LEVEL CHILLER 


temperature Tc via a dissipative thermal contact term of the 
form Sc® Be = 11) (2|® + h.c.. Here Be is the corresponding 

bath operator (cf Sec. IIB). Likewise, the transitions |1) 

|3) and |2) |3), with frequencies cu/, and tu* - a)c, couple 

weakly to heat baths at temperatures Tf, and T„, respectively. 
The temperatures are ordered asT„ > > Tc- 

One may compute the energy currents Qae{w,h,c] flowing 
from each of the heat baths into the system from its (non¬ 
equilibrium) steady state g(oo) [17]. Provided that 


OJe ^ ^c,rev 


CD/, 


Tc(Tw - T/,) 

Th(T„-Tey 


( 2 ) 


one has Q„ > 0, Q/, < 0 and Qc > 0. That is, net heat is 
extracted from the cold bath at temperature Tc and dumped 
into the hot bath at T/,, with the assistance of extra energy 
coming from the work bath at T„. This realizes a quantum 
absorption chiller [17, 18]. 

Unless cDc - Wc.rev the operation of the chiller is always 
irreversible, meaning that the total entropy production AS = 
-QwITw-Q/,IT/,-QcITc > 0 is strictly positive. This may be 
understood as a consequence of the ‘imperfect thermal con¬ 
tact’ between each transition and its heat bath. Indeed, let us 
define the internal temperatures Tq, [14, 36] as 


tD/i CDc (D/i CDc 

fw ~ / n \ ’ ~ / „ \ 5 fc — / n \ ’ 

log(g) log(g) log(g) 

where the p, stand for the steady-state populations in the en¬ 
ergy basis. The equality Tq, = only holds when cDc - Wc.rev 
[14]. In that limit AS - 0 and, as a result, the coefheient 
of performance e = Qc/Qw of the chiller saturates to the 
Carnot COP ec = Te(T„ - T/,)/[T„(Th - Tc)] [13, 14]. Other¬ 
wise, heat transfer is unavoidably irreversible due to the mis¬ 
match between Tq, and T^, and the COP evaluates instead to 
s = Tc(t„ - n)/[T„(T/, - Tc)] < Sc- 

Unlike the endoreversible three-level maser, the minimal 
four-level model described below does contain the relevant 
physical ingredients producing internal dissipation and heat 
leaks and therefore, it is more akin to a real absorption chiller. 


B. The four-level chiller and its steady-state solution 

Let us consider the following Hamiltonian 

//4 = +H[ = CDe m {b\ + \c) <c|) + CDh \d) {d\ 

+ g{\b){c\ + \c){b\). (4) 


A. Preliminaries: The three-level chiller 

Let us start by looking at the paradigmatic three-level ab¬ 
sorption chiller [17], with Hamiltonian 

//3 =Wc|2)<2|+a;,J3)<3|. (1) 

In all what follows h - ks - 1. The transition |1) <-» |2), 
of frequency cDc (< cu/,), is weakly coupled to a heat bath at 


We couple the transition \a} <-» \b) to the cold bath through 
a thermal contact operator of the form S c®Bc = \a) {b\ ®Bc + 
h.c, the transition \a) <-» \d), to the hot bath via a term like 
S h®Bk = \a) {d\ ®B/, + h.c., and finally |c) <-» \d), to the work 
bath with s |c) {c/| ®B„ H-h.c. (see Fig. 1). In all what 

follows, we assume g ^ cd/,. 

The eigenstates of H 4 are |1) s \a}, |2) s 2“'^^(|h) - |c)), 
|3) s 2“'^^(|h) -I- |c)), and |4) s \d}, and have energies 0, cDc - 
g, CDc + g, and CD/,, respectively. Henceforth, while the hot 
bath only acts on the transition |1) «-> |4), the work bath is 
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■\d) 

|c) 

\b) 

l«) 



Also since the dissipators X,a are in the LGKS form [40, 41], 
one can show that the irreversible entropy production is non¬ 
negative {AS > 0), as should be expected [42]. 

We now have at our disposal all the elements to carry out a 
full thermodynamic analysis of the four-level chiller. We shall 
use them below for its diagnosis and eventually, for devising 
a strategy to suppress the internal losses and heat leaks. 


FIG. 1: Sketch of the four-level chiller (left). The transitions cou¬ 
pled to each heat bath (work, hot and cold) are indicated by labelled 
arrows. The ^-dependent term in is depicted by a circle. On the 
right, the coupling scheme between the energy eigenstates is shown. 
The corresponding energies are indicated on the right-hand side. 


responsible for both the |2) <-> |4) and the |3) <-> |4) processes. 
Likewise, the cold bath acts on |1) <-> |2) and |1) <-> |3). 

We will assume bosonic reservoirs with Hamiltonian - 
Yjfi (a 6 (w, h, c]), where stands for the annihila¬ 

tion operator of mode in bath a. Hence, we may explicitly 
write Ba - g^iaa/j + where the coupling strengths are 
taken to be yield an Ohmic spectral den¬ 

sity [37]. This is the case for e.g. blackbody radiation. The 
parameter y sets the dissipation time scale td ~ 7 ^ ■ 

We also adopt the usual Born-Markov and rotating-wave 
approximations, implying that To is by far the largest time 
scale in the problem. Under this assumption one may read¬ 
ily derive an equation of motion for the state p(f) in the stan¬ 
dard Lindblad-Gorini-Kossakovski-Sudarshan (LGKS) form 
Q - {-Lw + -Ch + -Cc)q [38, 39], where the dissipators X.a are 
given by 

-Ca Q = ^a,a> ~ 2 7 (?}^j ’ (5) 

and (■, ■]+ stands for anti commutator. The non-Hermitian 
jump operators result from the decomposi¬ 

tion of Sa as eigenoperators of [37] (i.e. Sa = 
such that {Ha^Ac^u} - There is only one of them 

associated with the hot bath {A^ i^^^ - |1)(4|), which con¬ 
tributes to the corresponding dissipator with two terms, at fre¬ 
quencies ±LOh- On the other hand, the work bath adds four 
terms, at frequencies ±(a>„ + g), with the associated operators 
Aw.oj„.+g - |2)(4| and = |3)(4|. Here, we have de¬ 

fined cl)„ = (jJh - tOc- The same happens with the cold bath, 
that contributes with terms at frequencies +( 0 ;^ ±.g), for which 
A.,^,+^ = |l)<3|andA,,^^_, = |1)<2|. 

To ensure thermalization, the relaxation rates r„ must sat¬ 
isfy the detailed balance condition Ta,-ulTa^u - For 

blackbody radiation these are of the form = 7 W ^[1 -H 
riaico)], where tiaicxj) - [exp(cu/rQ.) - 1 - 1]“' is the bosonic oc¬ 
cupation number [37]. 

With all the above in mind, one only needs to impose sta- 
tionarity [i.e. {X.w + -Ch + -Cc)q{°°) = 0] to obtain the non¬ 
equilibrium steady state p(oo) of the four-level system. The 
heat currents can be defined as = trj// 4 XQ.p(oo)) [17]. 
Obviously, = 0 as a consequence of the stationarity. 


III. INTERNAL DISSIPATION AND HEAT LEAKS 

Heat devices are usually characterized by their power- 
versus-performance characteristic curves. Endoreversible 
models feature a distinctive open characteristic curve: When 
finite-rate heat transfer losses dominate, both the output power 
and the COP vanish. In the opposite (reversible) limit, the 
performance saturates to the Carnot COP, albeit at vanishing 
cooling load [see orange curves in Fig. 2(a)]. On the con¬ 
trary, more realistic devices affected as well by internal dis¬ 
sipation and heat leaks exhibit typical closed characteristic 
curves in which power and efficiency always vanish jointly 
[43]. The closed characteristic curves of our four-level ab¬ 
sorption chiller, depicted in Fig. 2(a) in blue, are thus a smok¬ 
ing gun evidence of additional irreversibility. Notice as well 
that increasing the parameter g seems to accentuate this effect. 
The aim of this section is to propose a meaningful breakdown 
of the total heat currents Q,y so as to identify the sources of the 
internal dissipation and parasitic heat transfer. 


A. Underlying stochastic dynamics 

Given a quantum master equation like the one above, one 
can always tackle the calculation of any physical quantity via 
the Monte-Carlo wave function (MCWF) method [44], that 
is, as an ensemble average over a sufficiently large number 
of stochastic trajectories. These are generated by a suitable 
piece-wise deterministic dynamics acting on pure states drawn 
from a given initial probability distribution. The MCWF 
method is not only technically convenient for computations 
involving large dissipative systems [37] but it may also pro¬ 
vide physical insights which would be hard to gasp from av¬ 
eraged quantities. 

In particular, the static viewpoint of the system frozen 
at q{oo) may be replaced by a dynamical one: that of 
a large number of time-evolving energy eigenstates drawn 
from the stationary distribution {pi,p2,pz,pA] = diagp(oo). 
Note that the state is written here in the energy eigenbasis 
(|1), |2), |3), |4)], rather than (|fl), \b ), |c), |^f)]. The dynamics 
along each trajectory would then consist of a sequence of in¬ 
dependent random jumps between the different energy eigen¬ 
states of the system, at the constant jump rates 

From the specific dissipative couplings of our four-level 
chiller, as depicted in Fig. 1(a), it is easy to see that any closed 
stochastic trajectory (i.e. starting and finishing at the same 
energy eigenstate), either falls into one of the following six 
categories, or is a permutation of combinations thereof: 
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FIG. 2: (a) Characteristic curves of the four-level chiller for g = jO. 1,0.3,0.5) (closed blue curves). The reservoir temperatures are T„ = 9, 
T;, = 8, Tc = 7, and the hot frequency is set to oi;, = 6 (in arbitrary units). To generate these curves, the cold frequency oj^ was increased 
until the device stopped cooling. Decoupling the device from all modes of the work bath at frequencies above e.g. w„, = oj/, - ojc completely 
suppresses the losses. The characteristic curves then become open like those of endoreversible devices (orange), and the COP saturates to 
its maximum possible value of sc = 7/9. (b) Individual contributions A5^ (solid), AS“ (dashed) and AS’'^“’‘ (dotted) to the total entropy 
production A5 = A5^ -l- A5“ AS*'^'*, in units of kg, around w^.rev for g = 0.1. The rest of parameters are the same as in panel (a). See 

discussion in Sec. Ill B 


(Cl) ID 
(C2) ID 
(C3) ID 
(C4) ID 
(C5) ID 
(C6) ID 


(i)c+8 i^w-g -COfi 

13) ^ |4) ^ ID. 

cold work hot 


CJh -(^w-g) -U^c+g) 

|4) -—24 |3) ^^4 ID. 

hot work cold 


(jJc-g (^w+g -liih 

|2)-D |4) ^ ID. 

cold work hot 


Lt)u -(a>„.+£) -iCJc-g) 

-4 |4) > |2) > ID. 

hot work cold 


(jJc-g (^w-^g ~{<^w~g) ~{iJ^c-^g) 

^ |2)-4 |4)-^ |3)-^ ID. 


cold 

{^c+g 

cold 


work 


work 


cold 


<j^w~g -(a>vv + e) -{Ci)c-g) 

D3) ^ |4) —-4^ |2) ID. 


work 


work 


cold 


The energy gained by the four-level system in each jump 
is indicated above the arrows and the corresponding baths are 
identified with the labels underneath them. 

Processes Cl and C3 are conventional thermodynamic 
cooling cycles with a net energy transfer from the cold to the 
hot bath [17], while C2 and C4 are their complementary heat¬ 
transforming cycles. On the other hand, the cycles C5 and 
C6 realize direct energy transfer between the work and the 
cold baths with no other effect on the hot reservoir, and are 
thus at the origin of the heat leaks. When it comes to pro¬ 
cesses C1/C2 and C3/C4 these may either work cooperatively 
or compete with each other. For g < ojc < coc.rev - g the rates 
of Cl and C3 are both larger than those of C2 and C4; that is 


Tc, -(oj^-hg) Tw, -(oj^i-g) h/i, ojf^ > Fc^ c>c+g hw, T/j, (ha) 

-(ajc-g) oih ^ oic-g oJw+g —ojh • (hb) 


Hence, the two processes cooperate producing net cooling, 
pretty much like the different stages in a multi-effect quantum 
absorption chiller [22]. However, in the range Wc.rev - g < 
aic < oJc^wv + g, the rate of C2 exceeds that of Cl while C3 
is still more likely to occur than C4. This means that the 
couple C3/C4 cools on average, as opposed to C1/C2, which 
produces net heating. This competition between cooling and 
heating mechanisms makes it impossible for the chiller to op¬ 
erate reversibly, and could be tagged ‘internal dissipation’ as 
it provides the leading contribution to the irreversible entropy 
production at ojc = a>c,rev (see Fig. 2(b) and the discussion in 
Sec. IIIB). For coc > a>c,iev +g, both C1/C2 and C3/C4 operate 
as heat transformers on average so that cooling is no longer 
possible. This is exactly the same mechanism that renders a 
periodically-driven three-level engine irreversible [10]. 

We shall formalize this intuition below, by suitably decom¬ 
posing the steady-state heat currents of our four-level chiller 
into contributions corresponding to two detuned elementary 
endoreversible stages, and an additional purely irreversible 
heat leak stage operating in parallel. 


B. Breakdown of the steady-state heat currents 


Let us introduce now the steady-state heat currents (34 s 
(Uy.. + g)I±, = -C0hl±, <34 = (coc + g)I±, and (3J,“'' = 

_(3ieak = g where the rates J±,ieak are given by 
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(7d) 








TlV, 



-r„ 


■+g ~ ^H-g > 



One can easily check that (3„ = for all 

three baths. The thermodynamic consistency of these newly- 
introduced stages contributing to the total steady-state heat 
currents, follows from the fact that they individually sum up 
to zero (i.e. - Yja^a - = 0)> and also sat¬ 

isfy the second-law-like inequalities AS^ = - YiaQal^a ^ 0, 
A5- = - > 0 and = - Ya^a^'^ITa > 0. 

We may associate the stage {(3,^) to the combined pro¬ 
cesses C1/C2, as the sign of 1+ is determined by the imbal¬ 
ance between the stationary rates of processes Cl and C2 [cf. 
Eqs. (7a) and (6a)]. Similarly, the sets of steady-state heat cur¬ 
rents {) and {) may be linked to processes C3/C4 and 
C5/C6, respectively. 

The ‘plus stage’ behaves as an absorption chiller (i.e. > 

0,Qt< 0 and > 0) within the cooling window 0 < oJc < 
Wc.rev - g, while the ‘minus stage’ does, provided that g < 
OJc < OJc,rev + g [21]. On the other hand, the heat leak rate 
Jieak is always negative since > Tc- This entails a parasitic 
heat current flowing directly from the work bath into the cold 
bath. Since we set g « w/„ the heat leaks will be typically 
negligible in comparison with the other contributions. Also 
note that as g increases, the cooling window for the plus stage 
shrinks while, for the minus stage, it is instead shifted as a 
whole towards larger frequencies. Thus, if we were to allow 
g to become larger than 0Jc,rev, the plus stage would not be 
capable of producing net cooling at any oJc, and the four-level 
chiller could ultimately lose its ability to cool. 

Just like in the case of the three-level chiller, the internal 
temperatures t^, and dehned from the steady-state pop¬ 
ulations {pi,p 3 ,p 4 ) in the energy basis as 




OJc 




( 8 ) 


coincide with the equilibrium temperatures of the reservoirs 
only at oJc = Wc.rev - g = tOccey At that point, the plus stage 
operates reversibly (AS^ - 0), i.e. at the Carnot COP. Sim¬ 
ilarly, dehning the internal temperatures from {pi,p 2 ,P 4 } 
one can see that the equalities and AS -Q hold at 

Wc = We, rev + ^ = tOcrev Howevef, the irreversible entropy 
production associated to the heat leak stage AS never van¬ 
ishes. That would require both r* = Tc and r* = T„ to hold 
simultaneously, which is not possible for any g 4 0. 


Fig. 2(b) shows the behaviour of the different shares of the 
total irreversible entropy production {AS^, AS^ and A5'®^) 
in the neighbourhood of oJc - 0 Jc,rev We can see that even if 
a non-vanishing alone is enough to render the device 

overall irreversible, the leading detrimental effect comes from 
the mismatch between and Also note that the 

larger the g, the larger the mismatch and, consequently, the 
bigger the resulting irreversibility at 0 Jc,rev This explains why 
increasing g lowers the maximum attainable COP in Fig. 2(a). 


C. Thermal short circuit: A different type of heat leaks 


Generalizing the above, we can say that a thermodynamic 
cooling cycle will feature irreversibility provided that two (or 
more) energy-basis transitions at different frequencies are al¬ 
lowed by the coupling to two (or more) heat baths. The dis¬ 
sipative transitions of such system may always be grouped to 
form coupled detuned three-level stages, which gives rise to 
both internal dissipation and heat leaks as illustrated above. 

An essentially different scenario is discussed in Ref [17], 
where irreversibility was modeled as the parasitic coupling 
of e.g. the cold bath to the work transition in a three-level 
maser. This would provide the chiller with an alternative dis¬ 
sipative path besides the conventional cooling and heating cy- 

cold work hot 

cles |1) «—» |2) «-» |3) «—> |1); that is, the shortcut 

work cold 

|2) <-> |3) <—> |2). As a result, heat could leak between 

the work and cold bath without producing any cooling. While 
the former effect is mere a by-product of internal dissipation, 
rooted in the interactions facilitating energy transfer within 
the working substance, the origin of this latter type of heat 
leaks, which we may tag ‘thermal short circuit’, lies simply 
in the lack of proper insulation at the thermal contacts of the 
device. Realizable nanoscale absorption chillers will typically 
be multi-terminal interacting devices, thus suffering from in¬ 
ternal dissipation and heat leaks. However, thermal short cir¬ 
cuit effects may also be unavoidable in practice, and could 
play a major role in the buildup of irreversible entropy pro¬ 
duction. 
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FIG. 3: (a) Sketch of the three-qubit absorption chiller and its tran¬ 
sition diagram. The eigen-energies are shown on the right, while the 
kets on the left denote the energy eigenstates following the notation 
of Ref. [21]. The transitions driven by the work, hot and cold baths 
are indicated with red, orange and blue arrows, respectively (see text 
for details), (b and c) Sketch of the two alternative configurations 
of the four-level chiller with their corresponding transition diagrams 
(see text for details). 


D. Suppression of internal dissipation and heat leaks 

From the above analysis it follows that the overall opera¬ 
tion of the four-level chiller can be made endoreversible by 
inhibiting either the plus or the minus stage. Let us imagine 
for instance that the spectrum of the work bath could be tai¬ 
lored into a Debye shape, with a high frequency cutoff 
anywhere between io„ — g and aj„ ■+ g. The corresponding 
excitation-relaxation rates would then read [1 + 

n„(a>)] 0(w„,_inax - |w|), where ©(■ ■ ■) stands for the Heavi¬ 
side step function. As a result, the stationary rates associated 
with the transition |2) <-4 |4) would vanish, thus precluding 
processes C3/C4 and C5/C6, and cancelling their contribu¬ 
tions to the total irreversible entropy production. Such engi¬ 
neered four-level chiller can operate reversibly at a>c = rev’ 
as showcased by its open characteristic curve, plotted in or¬ 
ange in Fig. 2(a). 

It is important to stress that the Markov approximation, 
which is central to our analysis, assumes that the decay of 
the two-time correlations of the baths is sufficiently fast com¬ 
pared with the dissipation time-scale. This assumption is jus- 


tihed so long as the spectral density remains approximately 
flat in the neighbourhood of the relevant Bohr frequencies 
of the system. Thus, for small g, the cutoff would be typ¬ 
ically too close to lJh, - g for the Markov approximation to 
be valid. As a result, both the actual time-evolution and the 
subsequent stationary state might differ significantly from the 
ones predicted by the LGKS quantum master equation. A rig¬ 
orous treatment of such problem requires more sophisticated 
tools, like higher order perturbative quantum master equa¬ 
tions [45, 46], non-Markovian stochastic Schrodinger equa¬ 
tions [47-49], hierarchical Heisenberg equations [50, 51], or 
other non-perturbative methods [52-54]. 

While the quantitative assessment of the performance of 
irreversible devices operating between engineered reservoirs 
is definitely worthy of investigation, it lies beyond the scope 
of this paper. Nevertheless, from the simple physical picture 
sketched above one may still propose practical means to min¬ 
imize the detrimental effects of heat leaks and internal dis¬ 
sipation. For instance, rather than connecting the transition 
|c) <-^ \d) directly to the work bath, one could interpose a har¬ 
monic oscillator of frequency atf between the two, so as to 
realize a frequency filter; The spectral density would then ef¬ 
fectively look like a skewed Lorentzian around tJf [55]. This 
reservoir-engineering technique has been recently considered 
in the context of nanoscale heat devices to produce spectrally- 
separated heat baths [26, 56]. Adjusting the detuning u„-a)f 
and the transition-oscillator coupling, one can easily find in¬ 
stances in which ^ 1, which strongly re¬ 

duces the rate of the minus stage in favour of the plus stage, 
and may thus allow for signihcantly larger cooling COPs. 


IV. THE THREE-QUBIT CHILLER 

We now turn our attention to a more complex three-terminal 
absorption chiller [18, 57]. These models are known to op¬ 
erate irreversibly as they feature closed characteristic curves 
[10, 20]. In particular, we consider three interacting two-level 
atoms, each one coupled to a heat bath, again at temperatures 
T„ > Th > Tc- As already mentioned, there exist experi¬ 
mental proposals for the implementation of such systems on 
superconducting qubits [31] or quantum dots [30]. 

More specifically, the model Hamiltonian reads 

//g — Utn;|lv,))(lvy| -H OJfi I 1/j) ( l/ll + (jJc I Ic) { Icl 

+ ^(|lw0/,lc)(0„U0c| H-h.C.), (9) 

in the eight-dimensional three-qubit Hilbert state space ® 
'K/, ® "Ke, where |0 q.) stands for the ground state and |1 q,), for 
the excited state [57] [see Fig. 3(a)]. The work frequency is 
set by the resonance condition a)„ = ojf, - a)c- The system- 
bath couplings may be written as Yja ® where 
is the Pauli x-matrix in 'Ha- Details on the derivation of the 
LGKS quantum master equation for this model can be found 
in Ref. [20]. 

It is crucial to remember that even if the terminals are lo¬ 
cally addressed by independent baths, the resulting dissipa- 
tors are global due to the three-body interaction term. It is 
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precisely the ensuing ‘delocalized dissipation’ that ultimately 
makes the device irreversible [20]. This hints that heat leaks 
could be at least partly responsible for the irreversibility. As 
pointed out in Ref. [58], an incorrect local modelling of the 
dissipation would not only fail to capture relevant physics of 
the problem, but might even predict physically unacceptable 
solutions which violate the second law of thermodynamics. 

A. Breakdown of the three-qubit chiller 

In Fig. 3(a) all the dissipative transitions between the en¬ 
ergy eigenstates of the three-qubit chiller are depicted with 
coloured arrows. Red stands for transitions coupled to the 
work bath, orange for those coupled to the hot bath, and blue 
for those of the cold bath. This device features nine open de¬ 
cay channels allowing for a total of eighteen dissipative tran¬ 
sitions between the eight energy eigenstates. Even if the dis¬ 
sipative dynamics of such model may seem quite involved at 
first glance, we will illustrate in what follows how our quali¬ 
tative understanding of the four-level chiller can greatly help 
in the system diagnosis and the subsequent performance opti¬ 
mization. 


1. Variations of the four-level chiller 

To proceed further, we must discuss first the two possible 
variations of the irreversible four-level chiller. The first one is 
illustrated in Fig. 3(b). Its Hamiltonian reads 

//; = w, \b} {b\ + cohilc) <c| + \d} {d\) + g(|c> {d\ + \d} <c|), (10) 

while the system-bath coupling operators are 5 c = \a) {b\ -H 
\b){a\, Sh - \a){d\ -H \d){a\, and S„ - |/7)<c| H- \c){b\. Re¬ 
peating the analysis from Sec. Ill for this alternative model, 
one finds that the competition between two detuned endore- 
versible stages is still responsible for most of the irreversible 
entropy production, the main difference being that the para¬ 
sitic heat current leaks from the work bath into the hot bath, 
rather than into the cold bath. On the other hand, the model 
depicted in Fig. 3(c) corresponds to the Hamiltonian 

= to, \c) {c\ + oJh \d) {d\ + g(\a} {b\ + \b} <fl|), (11) 

and system-bath coupling operators Sc = \b) (c| -t- |c) {b\, Sh - 
\a) {d\ H- \d) {a\, and = |c> {d\ -H \d) (c|. In this configuration 
the chiller features heat leaks from the hot bath into the cold 
bath instead. 


2. Diagnosis of the three-qubit chiller 

Generically, the dissipative transitions of any absorption 
chiller may be grouped into (possibly detuned) coupled three- 
level stages. We know from Sec. Ill B that the mismatch in the 
cooling windows of these stages gives rise to internal dissipa¬ 
tion, and that the frequency difference between pairs of transi¬ 
tions coupled to different baths facilitates heat leaks. One may 


easily identify the direction of the leaks by just pairing off 
these elementary three-level stages to yield either of the three 
configuration of the irreversible four-level chiller. In particu¬ 
lar, notice that in Fig. 3(a) we have split the dissipative tran¬ 
sitions of the three-qubit chiller into six coupled irreversible 
cycles, labelled with lowercase Roman numerals. One recog¬ 
nizes the transition diagram of Fig. 1 in stages (i) and (ii). In 
turn, stages (iii) and (iv) correspond to the configuration of 
Fig. 3(b), and stages (v) and (vi), to that of Fig. 3(c). Hence, 
we can conclude that this device suffers from all three types 
of pairwise heat leaks. However, for small g, it is again the 
mismatch in the cooling windows of the detuned three-level 
stages which contributes with the largest share to the total ir¬ 
reversible entropy production. 

Performance optimization would demand, in this case, to 
engineer the coupling to both the work and the hot baths. For 
instance, setting high-frequency cutoffs anywhere between 
(i)h - g and coh in the spectrum of the hot bath, and between 
Ww - g and (jOk in the spectrum of the work bath, suppresses 
all three-level stages except for the two that connect states [|3), 
|5), |8)], and (| 1), |4), |7)]. These are resonant and uncoupled, 
and the resulting device is therefore strictly endoreversible. 

We have thus showcased how one may diagnose a complex 
autonomous heat device, and even propose generic strategies 
to optimize its performance, without ever having to solve for 
its steady-state heat currents. 


B. COP at maximum power of irreversible devices 

The ultimate efficiency at maximum power that can be 
achieved by nanoscale heat devices has been object of intense 
study [19, 59-62]. In particular, for endoreversible single and 
multi-stage devices it has been established that the COP at 
maximum cooling load (e,) is tightly bounded from above as 
a result of the constraints imposed by the specific system-bath 
interaction mechanism [22, 23]. In particular, for unstructured 
bosonic baths, this bound is only set by the dimensionality 
of the heat reservoirs (e.g. e* < |ec for 3D baths) [21]. 
Whether the performance bound also applies to irreversible 
(‘non-ideal’) chiller models such as the three-qubit chiller was 
investigated via global numerical optimization in Ref. [20]. 
We are now in the position to give simple arguments to justify 
why the bound overestimates the maximum performance for 
any finite internal interaction strength. 

For simplicity, we turn our attention to the four-level chiller, 
though similar arguments apply as well to the three-qubit 
model. As pointed out in Sec. IIIB, the contribution of the 
heat leak component to the total steady-state heat currents and 
irreversible entropy production is rather small, so that one may 
approximate (3 q. ^ (3^ -H (3^;. In particular, the power of each 
stage Qf will peak at a different frequency and feature 
a COP of et = < |ec- However, the to¬ 

tal cooling load <3c will peak at some u)c,* + tof^, so that one 
has the strict inequalities e*(wc,*) < |sc- The overall COP at 


maximum cooling power may be thus approximated as 


e, - ^-e+(We_*) + ^-e_(Wc,*) < -rSc, (12) 


*) 




which is strictly below the endoreversible bound for unstruc¬ 
tured baths. In fact, is it easy to check that increasing the 
interaction strength g gradually lowers the actual bound on 
e*. Of course, if one is allowed to optimize e* globally (i.e. 
over all parameters), the endoreversible bound will tighten for 
vanishingly small interactions g. 

This is a neat example of why one should be careful not 
to extend results obtained within the endoreversible approxi¬ 
mation to generic irreversible devices. Especial care is in or¬ 
der when dealing with performance optimization, as endore- 
versibility typically leads to gross overestimations of the COP. 


V. CONCLUSIONS 

In this paper we have identified and characterized the 
mechanisms producing internal dissipation and heat leaks in 
nanoscale autonomous energy-conversion devices. These ef¬ 
fects are severely detrimental to their energy-efficient oper¬ 
ation and relate closely to each other. They arise whenever 
it becomes possible for one or more heat baths to excite two 
or more transitions of the working substance, at different fre¬ 
quencies. To that end, we have introduced and solved a mini¬ 
mal model of four-level irreversible absorption chiller. In par¬ 
ticular, we have split its steady-state heat currents into three 
contributions: two of them which correspond to coupled en¬ 
doreversible three-level stages operating in parallel, and an¬ 
other one which facilitates direct heat leaks from one heat 
bath into another. We have shown that, since the two endore¬ 
versible stages are detuned, there is a range of parameters for 
which they compete with each other instead of working co¬ 


operatively. It is precisely the mismatch in the cooling win¬ 
dows of these two stages which ultimately produces the largest 
share to the total irreversible entropy production as the coeffi¬ 
cient of performance grows. We have generically termed this 
effect ‘internal dissipation’ to differentiate it from the mere 
bath-to-bath parasitic heat leaks. We have also proposed the 
suppression of either of the two constituent three-level cy¬ 
cles of this irreversible model, via reservoir-engineering tech¬ 
niques, as a means to increase its maximum achievable coef¬ 
ficient of performance. Finally, we have shown how, aided by 
this qualitative understanding about the origin of irreversibil¬ 
ity in autonomous heat devices, one may be able to classify 
the sources of irreversible entropy production of more com¬ 
plex models of absorption chiller, e.g. determine the direc¬ 
tion of the heat leaks; and even to develop an intuition about 
potentially successful reservoir engineering techniques, to be 
implemented so as to allow for larger performance. 

A quantitative measure of the actual improvement in chiller 
performance achieved, for instance, by introducing hne-tuned 
frequency hlters at the thermal contacts, calls for a more care¬ 
ful modelling of the dissipative interactions beyond the com¬ 
mon Markovian approximation. It is certainly in order to 
pursue this research direction further, as it could ultimately 
pave the way towards a theory of thermal engineering at the 
nanoscale, with foreseeable widespread applications to quan¬ 
tum technologies. 
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